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Variational study of the ^He hypernucleus is presented using a realistic Hamiltonian and a fully 
correlated wave function including AiV space-exchange correlations. Behaviour of A-separation 
energy (Ba) with two- and three- baryon potential strengths is thoroughly investigated. Solutions 
for these potential strengths giving experimental Ba are presented. 
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Recently, a variational study [l| of the A He hypernu- 
cleus has been performed with a realistic Hamiltonian 
and a fully correlated wave function (WF). The WF takes 
into account all relavant dynamical correlations induced 
by the two- and three- baryon potentials and the AiV 
space-exchange correlation (SEC) that arises due to the 
space-exchange potential. The findings of the investi- 
gation suggest that no realistic study ignoring SEC is 
fair as it significantly affects every physical observable 
like energy breakdown, A-separation energy (Ba), nu- 
clear core polarization (NCP), point proton radius and 
density profiles. The effect is found more evident in the 
AA He double- A hypernucleus Q. The ground-state en- 
ergy of the hypernucleus (E) or the A-separation energy 
(-Ba = -E*Ho — E) depends on the strengths of the poten- 
tials involved in the Hamiltonian. 

A realistic Hamiltonian H of the hypernucleus, is writ- 
ten as a sum of the Hamiltonians due to the nuclear 
core (NC) of the hypernucleus (H^c) and due to the 
A-baryon (Ha), 



H = 
Hnc — 

H A = 



Hnc + Ha, 

TnC + ^2 Vij + ^2 V ijk, 
i<j i<j<k 



(i) 

(2) 
(3) 



Here, subscripts i,j and k refer to nucleons. For the 
5 = sector, we use Argonne vjr NN potential [3[ and 
Urbana type NNN potential [J [j| , which successfully 
explain the nuclear energy spectra and are well estab- 
lished. However, for the S — —1 sector, AN and ANN 
potential strengths are yet to be determined. The depen- 
dence of energy on these strengths is the theme of this 
work. 



Variation in any of the potential strengths would di- 
rectly affect the expectation value of the respective po- 
tential, ft would also affect the WF as correlations are a 
solutions of these potentials. Moreover, there are sensi- 
tivities to operators among various terms of the Hamil- 
tonian and the correlation functions. Correlations like 
SEC bring-in changes in the density profiles which af- 
fect even the central pieces of the energy breakdown. The 
basic ingredients like these strengths, therefore affect the 
ground-state energy collectively. It would not be possi- 
ble to perform a proper study for a particular potential 
strength ignoring others. Thus, in order to pin down 
these potential strengths, we have to handle them all to- 
gether. The energy of A Hc is thoroughly investigated 
along these lines. Such a study of all the s-shell single- 
and double-A hypernuclei may help to pin down these 
strengths, which, in turn, may resolve the outstanding 
A = 5 anomaly (g, [7|, lS| with no additional effort. Thus 
the present investigation is a step forward to answer, 
(i) whether we can successfully reproduce the hypernu- 
clear energy spectra using these potentials without in- 
voking the underlying Quantum Chromodynamics? (ii) 
and whether there is a possibility of physical existence 
of a bound (/ = 0, J = 1 + ) AA H hypernucleus which 
has been recently conjectured This would be helpful 
in studying the heavy hypernuclei specially 20 A Pb whose 
core density matches the nuclear matter density. Hence, 
this would lead us to investigate, in detail, the physics of 
charmed and bottom hypernuclei [HI El, E3] as well as 
to the nature and structure of neutron stars. 

The ANN potential arises from projecting out E, A, 
etc., degrees of freedom from a coupled channel formal- 
This is written as a sum of two terms, Va%j — 
f- Vfcj, as in Fig. [Tj The dispersive potential 
arising from the suppression mechanism owing to 
HiV coupling [H 0, [HI. flrl. is written including 
the explicit spin dependence as [17| 
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V& = W D TZ(m„rA l )Tl(m^A j )[l + <r A • + a^/Q}. 

(4) 
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However, spin term is too weak for spin zero core nu- 
cleus. Here, W D is the strength. The V^y is a two-pion 
exchange attractive potential. Neglecting higher partial 
waves it is written as a sum of two terms representing p- 
and s-wave ir — A scatterings, — V Ai 
Ref. [HI . The explicit form of these potentials are 



VXij as in 



V Aij 



(5) 



and 



with 



C s ' Z(m- K r i A)Z(m n rj A )cri ■ f iA <Tj ■ r jA T, t ■ Tj 

c s oi %] 



(6) 
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FIG. 1: Diagarm representing V® NN and V t 
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Xm = (cta ■ cr i )F 7r (m w rAj) + 5AiT' 7r (m 7r rAi) (7) 

and 

Z(x) = fa{x)-T«(x)]. (8) 

It may be expressed as generalised tensor-tau type opera- 
tors (<T i -rAi)(o- i -rAj), (cr i -r Al )(cr J TA l ), (cri-r Aj )(crj-rAj), 
and {(Ti-o-j) followed by (t^ -Tj), thus has a strong tensor 
dependence. In the above expressions, S A i is the tensor 
operator, Y n (x) is the Yukawa function 

Y«{x) = — (9) 
a; 

and T w (x) is the one-pion exchange tensor potential 

ZV(aO = f 1 + - + 4) — ^(r). (10) 
\ x x J x 

Here, £y(r) and ^r(r) are short-range cut-off functions, 

Mr) = 4 /2 W=(l-e- cr2 ). (11) 

Here, c = 2.0 fm~ 2 is a cut-off parameter and subscripts 
i,j and A refer to two nucleons and a A in the triplet 
(Aij). The C p and the C s are the strengths of V Ai j and 
V^ij , respectively. The latter is a very weak term com- 
pared to the former. Its strength is not known experi- 
mentally. However, we may make a qualitative theoreti- 
cal estimate for it by comparing the strengths of modern 
NNN potentials obtained using SU(3) s ymm etry. Us- 
ing chiral perturbation theory, Friar et al. [191 ] have com- 
pared the modern AAA potentials, namely: (i) Tucson- 
Melbourn [20j], (ii) Brazil [nj|, (hi) Ruhr [|| and (iv) 
Texas [23| containing a c-term for n—N scattering. They 
also consider the Fujita-Miyazawa force [24| dropping s- 
wave pions and Urbana-Argonne model [5| with addi- 
tional isospin- and spin- independent components added 
to the Fujita-Miyazawa force. 

In the Tucson-Melbourne (TM) model, the s-wave 
NNN force is written in the form 

B(Tij,Tj k ){Ti-Tj,Tj-Tk}{(Sij+<ri-<rj)(Sjk + tTj-tTk)}, 

(12) 



where, B{r i3 -, r^) has several terms given in Ref. p|. Re- 
cently, this has been expressed retaining only the term 
with pion-exchange-range functions as [25| 

A*=(^) a a'< (13) 
°ijk = J~] Z(m n rij)Z(mi r rjk)o-i-Tijcr k -r kj Ti-Tk. (14) 

eye 

The parameter a', whose value ranges from — 0.51/m w 
to -1.87/77V, is listed in Ref. [ig[. The TM value a' = 
— 1.03/tott gives the strength \A S \ «0.8 MeV. However, 
in many others it is assumed to have a value of 1.0 MeV. 

Comparing TM model with the Eq. [5] for ANN poten- 
tial, one may write an identical structure for both s-wave 
AAA and AAA potentials as following 

C s O s Ajk = A s O ijk . (15) 

This directly relates C s in the strange sector to A s in 
the non-strange sector. Since A — N mass difference is 
small compared to the A — A mass difference, AAA po- 
tential of S = — 1 sector is str ong er than its non-strange 
counterpart AAA potential [ijf of S = sector. This 
provides stronger strengths in the case of ANN potential 
compared to the NNN potential. We, therefore, expect 
that the value of C s would be more than 1.0 MeV, and 
is taken to be 1.5 MeV. 

The charge symmetric AN potential [2(| H3| reads as 

v\N(r) = v (r)(l - e + eP x ) + (7; CT /4)T 2 (m 7r r)cr A ■ ct n . 

The hrst term includes direct potential (vo(r) — V c (r) — 
V2-k{t)) and space-exchange potential [evq (r) (P x — 1)). 
Here, e determines the odd-state potential, which is 
the strength of the space-exchange potential relative to 
the direct potential. Its estimate from the Ap forward- 
backward asymmetry is poor that ranges from 0.1 to 
0.38 The potential v c (r) = W c /[l+exp{(r-R)/ar}} 
is the Saxon- Woods repulsive potential, with W c = 2137 
MeV, R = 0.5 fm and a = 0.2 fm, and v 27 r = vT^(m^r) 
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TABLE I: AN potential strengths in units of MeV. 





Vs 


v t 


v = (v a +3ut)/4 


V a = V s — Vt 


vl 


6.33 


6.09 


6.15 


0.24 


112 


6.28 


6.04 


6.10 


0.24 


v3 


6.23 


5.99 


6.05 


0.24 



is the two-pion attractive potential. The constants, 
v = (v s + 3^t)/4 and v a = v s — v t , are respectively the 
spin-average and spin-dependent strengths, with v s ( t ) the 
singlet (triplet) state potential depth. 

We perform variational Monte Carlo study to calculate 
the ground-state energy, E = where * 

is the WF of the hypernucleus. The computational de- 
tails are available in Ref. p). For the spin-zero core nu- 
cleus, the expectation value of the s pin p art of the AN po- 
tential is neglig ibly small 0, H, S3 H| . So is the s-wave 
part of ANN potential. Moreover, correlations induced 
by them are too weak to offer any significant change in 
the energy. Thus, we choose a reasonable strength for 
these two. The energy is very sensitive to changes in the 
AN potential strengths v and e. They implicitly appear 
in the WF through the AN central and the SEC correla- 
tions. The ANN potential and its correlations involving 
C p and W D play an important role. Therefore, E or B\ 
are sensitive to the strengths v, e, C p and W D . 

The value of v w 6.15(5) MeV is found consistent with 
the low energy Ap scattering data [13] • We use three dif- 
ferent sets of v s and Vt , which give three different values 
of v and a constant v a as in Table HI referred to as vl, 
V2 and T73. For all these, we choose three values of s in 
the range from 0.1 to 0.38 as mentioned before. These 
are 0.1, 0.2 and 0.3. Results for these values are given in 
Tables HH HI] and HVl 

The correlations induced by different components of 
ANN potential is written using scaled pair distances (r) 
and a variational parameter 5 m as in Ref. [28[ , 



E 

m=D,P,S 



S m V m (r Ml r l3l r ]A ). (17) 



dV^/dW D 



The (Vj^j) obeys a linear behaviour, 
<9£YcW D =constant at a fixed e as it is not sensitive to 
the operators but only to SEC and hence to e. As is 
obvious from Eq. [T71 a change in the strength C p offsets 
the WF. We observe that along with its own correlation 
parameter S p a couple of other parametres are found to 
change with C p . But the repulsive ANN correlation, 
U® NN , remains invariant. We perform calculations for a 
wide range of C p starting from 0.5 MeV and upto 2.0 
MeV. Therefore, enhancements in the attraction due to 
increasing C p needs to be balanced by an appropriate 
increase in W D . For every independent calculation, we 
tune the WF afresh and adjust the repulsive strength 
W D in order to reproduce the experimental B\. For the 
new W D (new), we may easily obtain new S D (new) using 

8 D {old)W D {old) = 5 D (new)W D (new), (18) 
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FIG. 2: Curves show the set of strengths giving B e ^ p . The 
dotted, dashed and long dashed lines represent e=0.1, 0.2 and 
0.3, and the red, black and blue colors represent vl, v2 and 
T>3, respectively. 



as U® NN is constant with C p . Thus S D decreases in the 
same proportion W D increases. The same is not true 
in the case of S p , which is found to be alomst constant 
even if we multiply C p by a factor of 4. This is because 
of the sensitivity of V p i3 to its correlation, which is so 
strong that the attraction, (V^), increases more than 
12 times for the corresponding 4 times increase in C p . 
This quadratic behaviour is observed for all the e and 
all the v. The respective increase in the NC part of 
the energy (Enc) is about 4 MeV. The T\ and v\i also 
exhibit considerable change due to the variation in C p . 

Solutions of all these strengths reproducing B e ^ p are 
plotted in Fig. [5] Thus, following the curves one may 
find potential strengths that reproduce B\ '. For every 
v and at a fix value of C p , we observe a linear realtion- 
ship, dW D /de « c, which is a consequence of two other 
linear relationships: (i) dE/de = —dB^/de w ci (ii) 
and dVfij/dW = -dE/dW D = dB A /dW D ps c 2 . The 
slope, c = c\/c2 is found to increases with C p , but only 
slightly (Table |V|) . Curves representing different v are 
found to get closer at higher C p . Because, the depen- 
dence of energy on C p as well as on W D varies with 'v. 
To match an increase in the attractive C p we require a 
larger increase in the repulsive W D for smaller v. 

This study has established the value of parameters 
appearing in the two- and three-body potentials in the 
strange sector. Furthermore, range of variations of these 
parameters is now, atleast roughly, known. In order to 
establish these parameters over a range of hypernuclci, 
numerous light and heavy hypernuclei have to be stud- 
ied. This would help us to understand the variation of 
AA*" and ANN potentials as a function of density. Such a 
study would be in parallel to that of NN and NNN po- 
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tentials in nuclear systems where density variation plays 
an important role. A proper understanding of hypernu- 
clei would help us to clarify the behaviour of the nuclear 
forces in the strange sector. Ultimately, a detailed study 
may lead us to a clarification of the role of QCD in de- 
termining the potential strengths. New results expected 



from Japan Hadron Facility would help to sort out these 
questions in the near future. 
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32/99 awarded to AAU by the Department of Science and 
Technology Government of India. FCK thanks NSERCC 
for financial support. 



[1] A. A. Usmani, Phys. Rev. C73, 011302(R) (2006). 

[2] A. A. Usmani and Z. Hasan, submitted to Phys. Rev. C. 

[3] R. B. Wiringa, V. G. J. Stoks, and R. Schiavilla, Phys. 

Rev. C51, 38 (1995). 
[4] B. S. Pudliner, V. R. Pandharipande, J. Carlson and R. 

B. Wiringa, Phys. Rev. Lett. 74, 4396 (1995). 
[5] J. Carlson, V. R. Pandharipande and R. B. Wiringa, 

Nucl. Phys. A401, 59 (1983). 
[6] R. H. Dalitz, R. C. Herndon and Y. C. Tang, Nucl. Phys. 

B47, 109 (1972). 
[7] A. Gal, Adv. Nucl. Phys. 8, 1 (1975). 
[8] E. V. Hungerford and L. C. Biedenhorn, Phys. Lett. 

142B, 232 (1984). 
[9] K. Ahn et ai, Phys. Rev. Lett 87, 132504 (2001). 
[10] A. A. Tyapkin, Sov. J. Nucl. Phys. 22, 89 (1976). 
[11] C. B. Dover and S. H. Kahana, Phys. Rev. Lett. 39, 1506 

(1977). 

[12] K. Tsushima and F. C. Khanna, Phys. Rev. C67, 015211 
(2003). 

[13] A. R. Bodmer, D. M. Rote and A. L. Mazza, Phys. Rev. 

C2, 1623 (1970). 
[14] A. R. Bodmer and D. M. Rote Nucl. Phys. A169, 1 

(1971). 

[15] J. Rozynek and J. Dabrowski, Phys. Rev. C20, 1612 
(1979); 23, 1706(1981); Y. Yamamoto and H. Bando, 
Suppl. Prog. Theor. Phys. Suppl. 81, 9(1985); Y. Ya- 
mamoto, Nucl. Phys. A450, 275c (1986). 

[16] Q. N. Usmani and A. R. Bodmer, Phys. Rev. C60, 
055215 (1999). 



[17] A. R. Bodmer and Q. N. Usmani, Nucl. Phys. A477, 621 
(1988). 

[18] R. K. Bhaduri, B. A. Loiseau and Y. Nogami, Anns. 

Phys. (N. Y) 44, 57 (1967). 
[19] J. L. Friar, D. Hiiber and U. van Kolck, Phys. Rev. C59, 

53 (1999). 

[20] S. A. Coon, M. D. Scadron, P. C. McNamee, B. R. Bar- 
rett, D. W. E. Blatt and B. H. J. McKellar, Nucl. Phys. 
A317, 242 (1979). 

[21] H. T. Coelho, T. K. Das and M. R. Robilotta, Phys. Rev. 
C28, 1812 (1983); M. R. Robilotta and H. T. Coelho, 
Nucl. Phys. A460, 645 (1986). 

[22] J. A. Eden and M. F. Gari, Phys. Rev. C53, 1510 (1996). 

[23] C. Ordonez and U. van Kolck, Phys. Lett. B 291, 459 
(1992); U. van Kolck, Phys. Rev. C49, 2932 (1994). 

[24] J.-I Fujita and H. Miyazawa, Prog. Theor. Phys. 17, 360 
(1957). 

[25] S. C. Pieper, V. R. Pandharipande, R. B. wiringa and J. 

Carlson, Phys. Rev. C64, 014001 (2001). 
[26] A. R. Bodmer and Q. N. Usmani and J. Carlson, Phys. 

Rev. C29, 684 (1984). 
[27] I. E. Lagaris and V. R. Pandharipande, Nucl. Phys. 

A359, 331 (1981). 
[28] A. A. Usmani and S. Murtaza, Phys. Rev. C68, 024001 

(2003). 

[29] A. A. Usmani, S. C. Pieper and Q. N. Usmani, Phys. 

Rev. C51, 2347 (1995). 
[30] A. A. Usmani, Phys. Rev. C52, 1773, (1995). 



TABLE II: Energy breakdown for the set vl (v = 6.15 and v a — 0.24). All quantities are in units of MeV except for e. 







e — 


0.1 






£ = 


0.2 






£ — 


0.3 






C p = .5 


C p = 1.0 


C p = 1.5 


C p = 2.0 


C p = .5 


C p = 1.0 


C p = 1.5 


C p = 2.0 


C p = .5 


C p = 1.0 


C p = 1.5 


C p = 2.0 


Ta 


8.56(3) 


8.97(3) 


9.83(3) 


10.46(4) 


8.23(3) 


8.76(3) 


9.38(3) 


10.02(4) 


7.94(3) 


8.39(3) 


9.01(3) 


9.58(3) 


v (r)(l - e) 


-16.08(6) 


-16.36(6) 


-16.78(6) 


-17.03(6) 


-13.69(5) 


-14.07(5) 


-14.37(5) 


-14.51(5) 


-11.55(4) 


-11.91(4) 


-12.06(4) 


-12.22(5) 


vo(r)eP x 


-1.56(1) 


-1.58(1) 


-1.61(1) 


-1.63(1) 


-2.97(1) 


-3.04(1) 


-3.09(1) 


-3.11(1) 


-4.27(2) 


-4.40(2) 


-4.43(2) 


-4.47(2) 


{\)v a Tl{r)<TA ■ (Ti 


0.014(0) 


0.015(0) 


0.017(0) 


0.017(0) 


0.012(0) 


0.013(0) 


0.014(0) 


0.014(0) 


0.009(0) 


0.010(0) 


0.011(0) 


0.011(0) 


VAi 


-17.63(6) 


-17.93(6) 


-18.38(6) 


-18.65(7) 


-16.65(6) 


-17.10(6) 


-17.42(6) 


-17.61(7) 


-15.80(6) 


-16.30(6) 


-16.48(6) 


-16.68(7) 




2.49(1) 


4.60(1) 


7.99(4) 


12.17(6) 


2.16(1) 


4.20(1) 


7.35(4) 


11.53(6) 


1.80(1) 


3.81(1) 


6.81(4) 


10.65(5) 


Vl ■■ 


-1.36(1) 


-4.40(2) 


-10.02(4) 


-15.99(6) 


-1.36(1) 


-4.54(2) 


-9.42(4) 


-15.35(6) 


-1.28(1) 


-4.16(2) 


-8.97(4) 


-14.38(6) 


s 3 

P S 


-0.025(0) 


0.017(1) 


0.066(1) 


0.100(1) 


-0.025(0) 


0.009(1) 


0.057(1) 


0.087(1) 


-0.030(1) 


0.009(1) 


0.050(1) 


0.077(1) 


^ Aij ^Aij ^Aij 


-1.38(1) 


-4.39(2) 


-9.95(4) 


-15.89(6) 


-1.39(1) 


-4.52(2) 


-9.36(4) 


-15.26(6) 


-1.31(1) 


-4.16(2) 


-8.92(4) 


-14.30(6) 


VAij = V Ai j + V A ij 


1.11(1) 


0.21(2) 


-1.96(2) 


-3.73(3) 


0.77(1) 


-0.33(2) 


-2.05(2) 


-3.73(3) 


0.49(1) 


-0.34(2) 


-2.11(2) 


-3.65(3) 


Va = VAi + VAi] 


-16.52(6) 


-17.72(6) 


-20.34(6) 


-22.38(7) 


-15.88(6) 


-17.43(6) 


-19.43(6) 


-21.34(7) 


-15.31(6) 


-16.64(6) 


-18.59(6) 


-20.33(7) 


Ea = Ta + Va 


-7.96(3) 


-8.75(4) 


-10.51(4) 


-11.92(5) 


-7.64(4) 


-8.67(4) 


-10.05(4) 


-11.32(5) 


-7.37(4) 


-8.26(4) 


-9.58(4) 


-10.75(5) 


Tnc 


117.59(15) 


118.42(15) 


118.69(15) 


119.26(15) 


116.97(15) 


117.84(15) 


118.42(15) 


118.84(15) 


117.51(15) 


118.43(15) 


118.07(15) 


118.53(15) 


vnn 


-134.65(14) 


-134.53(14) 


-133.18(14) 


-132.55(14) 


-134.36(14) 


-134.05(14) 


-133.41(14) 


-132.79(14) 


-134.99(14) 


-134.85(14) 


-133.42(14) 


-132.87(14) 


Vnnn 


-5.81(2) 


-5.98(2) 


-5.84(2) 


-5.66(2) 


-5.82(2) 


-5.97(2) 


-5.83(2) 


-5.67(2) 


-5.99(2) 


-6.15(2) 


-5.93(2) 


-5.76(2) 


Vnc = Vij + Vijk 


-140.46(14) 


-140.51(14) 


-139.02(14) 


-138.21(14) 


-140.16(14) 


-140.03(14) 


-139.24(14) 


-138.46(14) 


-140.98(14) 


-141.00(14) 


-139.35(14) 


-138.63(14) 


Enc = Tnc + Vnc 


-22.88(4) 


-22.10(5) 


-20.32(4) 


-18.95(5) 


-23.21(4) 


-22.18(5) 


-20.82(4) 


-19.52(5) 


-23.47(4) 


-22.58(5) 


-21.27(4) 


-20.11(5) 


E = E a + E nc 


-30.84(2) 


-30.85(2) 


-30.84(3) 


-30.87(4) 


-30.85(2) 


-30.85(2) 


-30.86(3) 


-30.84(4) 


-30.84(2) 


-30.84(2) 


-30.86(3) 


-30.86(4) 



Cn 



TABLE III: Energy breakdown for the set v2 (v — 6.10 and v a = 0.24). All quantities are in units of MeV except for e. 







£ = 


U.l 






£ = 


{J.Z 






£ = 


U.o 






C p = .5 


C p = 1.0 


C p = 1.5 


C p = 2.0 


C p = .5 


C p = 1.0 


C p = 1.5 


C p = 2.0 


C p = .5 


C p = 1.0 


C p = 1.5 


C p = 2.0 


Ta 


8.09(3) 


8.54(3) 


9.19(3) 


9.91(3) 


7.80(3) 


8.15(3) 


8.65(3) 


9.51(3) 


7.57(3) 


7.94(3) 


8.44(3) 


9.18(3) 


vo(r)(l - e) 


-14.61(5) 


-14.79(5) 


-15.13(5) 


-15.48(5) 


-12.54(5) 


-12.64(5) 


-12.94(5) 


-13.30(5) 


-10.59(4) 


-10.73(4) 


-11.04(5) 


-11.15(4) 


vo(r)eP x 

(iH2?(r)ff A • tTi 


-1.41(1) 


-1.42(1) 


-1.45(1) 


-1.47(1) 


-2.70(1) 


-2.72(1) 


-2.77(1) 


-2.83(1) 


-3.88(2) 


-3.93(2) 


-4.03(2) 


-4.03(2) 


0.007(0) 


0.007(0) 


0.008(0) 


0.009(0) 


0.007(0) 


0.008(0) 


0.008(0) 


0.009(0) 


0.005(0) 


0.005(0) 


0.006(0) 


0.007(0) 


VAi 


-16.02(6) 


-16.20(6) 


-16.57(6) 


-16.94(6) 


-15.23(6) 


-15.35(6) 


-15.70(6) 


-16.12(6) 


-14.46(6) 


-14.65(6) 


-15.04(6) 


-15.18(6) 


Vaij 


1.65(1) 


3.46(2) 


6.43(4) 


10.39(5) 


1.32(1) 


3.09(2) 


5.95(4) 


9.81(5) 


0.99(1) 


2.71(2) 


5.47(4) 


9.22(5) 


Vl- 


-1.33(1) 


-4.36(2) 


-9.00(4) 


-14.95(6) 


-1.23(1) 


-4.04(2) 


-8.08(3) 


-14.26(6) 


-1.19(1) 


-3.98(2) 


-7.61(3) 


-13.84(6) 


s J 

V Aij 


-0.028(0) 


0.014(1) 


0.049(1) 


0.087(0) 


-0.035(0) 


0.007(0) 


0.040(1) 


0.078(0) 


-0.032(0) 


0.009(0) 


0.036(1) 


0.074(0) 


V Aij ^Aij Vaij 


-1.35(1) 


-4.35(2) 


-8.59(4) 


-14.86(6) 


-1.27(1) 


-4.04(2) 


-8.04(4) 


-14.18(6) 


-1.22(1) 


-3.97(2) 


-7.80(4) 


-13.77(6) 


Vazj = V Ai j + V Ai j 


-0.30(1) 


-0.88(2) 


-2.53(2) 


-4.47(3) 


0.05(1) 


-0.94(2) 


-2.09(2) 


-4.37(3) 


0.49(1) 


-1.26(2) 


-2.33(2) 


-4.55(3) 


Va = VAi + VAi 3 


-15.72(6) 


-17.08(6) 


-19.10(6) 


-21.41(6) 


-15.18(6) 


-16.29(6) 


-17.79(6) 


-20.49(6) 


-14.69(6) 


-15.91(6) 


-17.37(6) 


-19.73(6) 


Ea = T a + V a 


-7.62(3) 


-8.54(4) 


-9.90(4) 


-11.51(5) 


-7.37(3) 


-8.14(4) 


-9.15(4) 


-10.98(5) 


-7.13(4) 


-7.97(4) 


-8.93(4) 


-10.54(5) 


Tnc 


116.82(15) 


116.53(15) 


117.36(15) 


118.09(15) 


116.73(15) 


116.68(15) 


117.19(15) 


117.55(15) 


116.81(15) 


117.03(15) 


117.16(15) 


117.40(15) 


vnn 


-134.22(14) 


-132.90(14) 


-132.59(14) 


-131.80(14) 


-134.40(14) 


-133.50(14) 


-133.12(14) 


-131.82(14) 


-134.43(14) 


-133.72(14) 


-133.12(14) 


-131.96(14) 


Vnnn 


-5.83(2) 


-5.92(2) 


-5.74(2) 


-5.62(2) 


-5.81(2) 


-5.90(2) 


-5.78(2) 


-5.60(2) 


-6.10(2) 


-6.21(2) 


-5.94(2) 


-5.73(2) 


Vnc = Vij + Vijk 


-140.05(14) 


-138.82(15) 


-139.31(15) 


-137.42(14) 


-140.21(14) 


-139.41(14) 


-138.90(15) 


-137.42(14) 


-140.53(14) 


-139.93(15) 


-133.12(15) 


-137.69(14) 


Enc = Tnc + Vnc 


-23.22(4) 


-22.31(5) 


-20.97(5) 


-19.32(5) 


-23.47(4) 


-22.71(5) 


-21.71(5) 


-19.87(5) 


-23.71(4) 


-22.89(4) 


-21.91(4) 


-20.29(5) 


E = E a + E nc 


-30.84(2) 


-30.85(2) 


-30.87(2) 


-30.84(3) 


-30.85(2) 


-30.85(2) 


-30.86(2) 


-30.85(3) 


-30.84(2) 


-30.87(2) 


-30.84(2) 


-30.84(3) 



C5 



TABLE IV: Energy breakdown for the set «3 (v = 6.05 and v a = 0.24). All quantities are in units of MeV except for e. 







£ = 


U.l 






£ = 


{J.Z 






£ = 


U.o 






C p = .5 


C p = 1.0 


C p = 1.5 


C p = 2.0 


C p = .5 


C p = 1.0 


C p = 1.5 


C p = 2.0 


C p = .5 


C p = 1.0 


C p = 1.5 


C p = 2.0 


Ta 


7.60(3) 


8.09(3) 


8.66(3) 


9.18(3) 


7.27(3) 


7.73(3) 


8.14(3) 


8.87(3) 


7.11(3) 


7.45(3) 


7.90(3) 


8.48(3) 


vo(r)(l - e) 


-13.26(5) 


-13.57(5) 


-13.69(5) 


-13.95(5) 


-11.20(4) 


-11.48(5) 


-11.64(5) 


-12.02(5) 


-9.58(4) 


-9.68(4) 


-9.90(4) 


-9.94(4) 


vo(r)eP x 

(iK2?(r)ff A - «Ti 


-1.27(1) 


-1.30(1) 


-1.30(1) 


-1.31(1) 


-2.40(1) 


-2.45(1) 


-2.48(1) 


-2.54(1) 


-3.49(2) 


-3.51(2) 


-3.56(2) 


-3.58(2) 


-0.003(0) 


-0.002(0) 


-0.002(0) 


-0.001(0) 


0.003(0) 


0.004(0) 


0.005(0) 


0.006(0) 


0.002(0) 


0.002(0) 


0.003(0) 


0.003(0) 


«Ai 


-14.53(6) 


-14.87(6) 


-15.00(6) 


-15.27(6) 


-13.60(6) 


-13.92(6) 


-14.10(6) 


-14.56(6) 


-13.06(6) 


-13.19(6) 


-13.48(6) 


-13.52(6) 


^Aij 


0.87(0) 


2.61(1) 


5.30(3) 


8.79(5) 


0.51(0) 


2.16(1) 


4.68(3) 


8.25(5) 


0.24(0) 


1.82(1) 


4.31(3) 


7.50(4) 




-1.20(1) 


-4.07(2) 


-8.53(4) 


-13.45(6) 


-1.11(1) 


-3.84(2) 


-7.49(4) 


-12.87(6) 


-1.06(1) 


-3.62(2) 


-7.24(4) 


-12.27(6) 


s 3 


-0.032(0) 


0.010(0) 


0.046(0) 


0.071(1) 


-0.037(0) 


0.003(0) 


0.030(0) 


0.062(1) 


-0.037(0) 


0.001(0) 


0.029(0) 


0.060(1) 


Vaij Vazj Vaij 


-1.23(1) 


-4.06(2) 


-6.48(4) 


-13.37(6) 


-1.15(1) 


-3.84(2) 


-7.46(4) 


-12.81(6) 


-1.09(1) 


-3.63(2) 


-7.21(4) 


-12.21(6) 


Vaij = Vajj + Vajj 


-0.37(1) 


-1.44(2) 


-3.19(2) 


-4.59(3) 


-0.64(1) 


-1.68(2) 


-2.78(2) 


-4.56(3) 


-0.85(1) 


-1.81(2) 


-2.90(2) 


-4.71(3) 


Va = «Ai + Vaij 


-14.89(7) 


-16.32(6) 


-18.19(6) 


-19.86(6) 


-14.24(7) 


-15.60(6) 


-16.89(6) 


-19.12(6) 


-13.91(7) 


-11.48(6) 


-16.38(6) 


-18.23(6) 


£a = T A + V A 


-7.29(3) 


-8.23(3) 


-9.52(4) 


-10.67(4) 


-6.97(3) 


-7.87(4) 


-8.75(4) 


-10.25(4) 


-6.80(3) 


-7.55(4) 


-8.48(4) 


-9.74(4) 


Tnc 


115.97(15) 


116.33(15) 


116.54(15) 


116.42(15) 


115.28(15) 


115.70(15) 


116.01(15) 


116.81(15) 


116.04(15) 


116.28(15) 


116.23(15) 


115.68(15) 


vnn 


-133.81(14) 


-133.07(14) 


-132.06(14) 


-131.98(14) 


-133.44(14) 


-132.81(14) 


-132.32(14) 


-131.73(14) 


-134.04(14) 


-133.41(14) 


-132.63(14) 


-131.08(14) 


Vnnn 


-5.73(2) 


-5.88(2) 


-5.80(2) 


-5.61(2) 


-5.73(2) 


-5.87(2) 


-5.75(2) 


-5.67(2) 


-6.03(2) 


-6.16(2) 


-5.95(2) 


-5.71(2) 


Vnc = Vij + Vijk 


-139.53(14) 


-138.95(14) 


-137.86(14) 


-136.59(14) 


-139.17(14) 


-138.68(14) 


-138.07(14) 


-137.40(14) 


-140.08(14) 


-139.57(14) 


-138.58(14) 


-136.79(14) 


Enc = Tnc + Vnc 


-23.57(4) 


-22.61(5) 


-21.32(4) 


-20.17(5) 


-23.89(4) 


-22.98(4) 


-22.09(4) 


-20.59(5) 


-24.04(4) 


-23.29(4) 


-22.36(4) 


-21.11(5) 


E = E a + E nc 


-30.85(2) 


-30.84(2) 


-30.84(2) 


-30.84(3) 


-30.85(2) 


-30.85(2) 


-30.84(2) 


-30.84(3) 


-30.84(2) 


-30.84(2) 


-30.85(2) 


-30.86(3) 



-J 



TABLE V: Variation of the slope, dW D /de, with C p and v. 



C p (MeV) 


TJl (MeV) 


v2 (MeV) 


v3 (MeV) 


0.5 


-0.016(1) 


-0.017(1) 


-0.017(1) 


1.0 


-0.017(1) 


-0.019(1) 


-0.019(1) 


1.5 


-0.019(1) 


-0.022(1) 


-0.023(1) 


2.0 


-0.021(1) 


-0.023(1) 


-0.024(1) 


2.5 


-0.022(1) 


-0.025(1) 


-0.026(1) 



